APPROXIMATION OF THE SEMIGROUP GENERATED BY 
THE ROBIN LAPLACIAN IN TERMS OF THE GAUSSIAN 

SEMIGROUP 
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Abstract. For smooth bounded open sets in euclidean space, we construct 
corresponding contractive linear extension operators for the space of continu- 
ous functions which preserve regularity of functions in the domain of the Robin 
Laplacian. We also prove a Trotter-like approximation for the semigroup gen- 
erated by the Laplacian subject to Robin boundary conditions in terms of 
these extension operators. The limiting case of Dirichlet boundary conditions 
is treated separately. 



1. Introduction 



Let il C be a smooth bounded open set. Here and in the following, "smooth" 
means "of class C°°", although the main results remain true under slightly milder 
regularity assumptions. On such a set we consider the (autonomous, homogeneous) 
diffusion equation 



Ut = Au, 

Oil 

-Q^{t,z) = -I3{z)u{t,z), 
u{0,x) = uo{x), 



on (0, oo) X Q, 

for t > and z G d^l, 
for a; e ri. 



(1) 



subject to Robin boundary conditions. Here uq e C{fl) is an arbitrary initial 
function, /3 is a non-negative smooth function on dfl which does not depend on 
t and ^ denotes the directional derivative of u along the outwards pointing unit 
normal of fl. We remark that this setting includes Neumann boundary conditions 
for (3 = 0. A mild solution of ([T]) is a function u £ C([0, oo); C(ri)) such that 
/p M(s)ds e D{Ar) and 

u{t) = Mo + Ar / u{s)ds 
Jo 

for every t > 0. Note that we have incorporated the Robin boundary conditions 

— + (3u = Qondn (2) 
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into the domain 



D{Ar) := \ ue H\Q)nc{n) 



/ \7u\7ip+ / Auip+ / uip /3da = ioT every if e H^{Q)\ , 
AffU :— Au 

of the Laplacian on Q subject to Robin boundary conditions. 

It is known that for every non-negative, bounded, measurable function /3 and 
every initial value uq G G{Q) there exists a unique mild solution to problem ([1]). In 
fact, Warma proved that under the above assumptions generates a Co-semigroup 
Tii{t) on C(il) [14], and it follows from the general theory of semigroups that 
then u{t) = Tii[t)u{) is the unique mild solution of the corresponding homogeneous 
abstract Cauchy problem 4, Proposition II. 6. 4]; note that Warma's proof remains 
valid for arbitrary non- negative functions j3 G L°°{d^l). 

If we want to calculate this solution numerically, a typical problem is how to 
handle the boundary conditions. To fix the general ideas, let H. — (0, 1) and assume 
that we want to apply an explicit finite difference method. Then one replaces 
the derivatives Ut and u^x by appropriate difference quotients and successively 
calculates approximations u{tn,Xj) of the exact solution u(t,x) by the relation 

u[tn+l,Xj) ~ u{tn,Xj) _ u{tn,Xj+i) - 2u{tn,Xj) + tt(tn,Xj-i) 
k ^ 

where tn — n ■ k and Xj — j ■ h for given small numbers k,h > 0. Note, however, 
that this cannot be directly applied to calculate u(tn+i,0) and u{tn+i, 1) because 
u(tn, —h) and u{tn, l+h) are not defined. For Dirichlet boundary conditions, we can 
assign u(t„, —h) :— u{tn, l + h) :— 0. On the other hand, for Neumann boundary 
conditions the situation is not as simple. One common technique is to use 

u{tn,~h) :— u{tn,h) and 7i(i„, 1 -|- /i) := u(t„, 1 — ft,) (3) 

in the calculations, which comes from a second order accurate approximation of the 
derivative at the boundary, see [TTl Section 8.3]. 

The aim of the article at hand is to justify the use of ([3]) for Neumann boundary 
conditions from a semigroup perspective, showing that the corresponding continu- 
ous method converges to the exact solution as /c ^ 0, and to extend it to the more 
general case of Robin boundary conditions. More precisely, we construct an exten- 
sion operator from C(r2) to Co(R^) which depends only on and /3 (but not on 
t) and resembles a continuous version of ([3]) if /3 = 0, such that Ep is a contraction 
and that EfjU is sufficiently regular whenever u g £>(Ar); we refer to Corollary [T2l 
for the precise statement. For operators Efs satisfying these two assumptions, we 
prove the Trotter- like (compare to [T2]) approximation result 

TR{t)u = lim {RGo{^)EaYu in 0(1!) for every u e C(II) (4) 

uniformly on [0,r] for every T > 0, where Go(t) denotes the Gaussian semigroup 
on Co(lR^) and R: Co(lR^) — > C(il) is the restriction operator Ru — u\n- This 
shows how Robin boundary conditions can be incorporated into a numerical solver 
such that the numerical solutions converge uniformly on [0, T] x fi, at least if error 
introduced by space discretization is neglected. 
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The article is organized as follows. In Section [2] we show how to represent a 
neighborhood of dfl in terms of the outwards pointing unit normal and recall some 
facts about the Laplacian. In Section [3] we construct the extension operator Ep 
related to the Robin Laplacian and prove its aforementioned two properties which 
ensure (j4]) as we show in Section |4l Section [5] deals with the limiting case /3 — > oo 
giving rise to Dirichlet boundary conditions. Finally, Section [6] summarizes the 
results. 



2. Notation and Preliminary Results 

It is well-known that for smooth boundary a neighborhood of dQ can be para- 
metrized by the outwards pointing unit normal v. Because certain features of the 
parametrization are needed later on, we state this result in the formulation we want 
to use and prove it for the sake of completeness. 

Proposition 1. Let S > and 

T:dQx {-S, (5) ^ R^, (p, t)^p + tu{p). 
For small 5 > 0, T is a smooth diffeomorphism onto a neighborhood of dO, in R^. 
Remark 2. 

(a) Let xq G dn be arbitrary. By definition of "smooth boundary" , fi can locally at 
xq be represented as the subgraph of a smooth function (p: [/ ^ R ([/ C R^~^) 
up to rotation. Assume for the moment that no rotation is needed. Then 



dnnv 



z 
viz] 



zeU} (5) 



for an open set V C R^. Thus z i-^ ( ip(^z) ) is a bijection of an open subset of 
R^-i onto a neighborhood of xq in dfl. Using these mappings for all x G dfl 
as charts we make dil into a smooth manifold built upon the subspace topology 
induced by R^. Thus we can look at T as a mapping from a manifold to a 
euclidean space. As usual we say that T is smooth if the composition T* of T 
with a chart is smooth, i.e., if 

is smooth as a mapping between euclidean spaces, 
(b) Using charts, the outwards pointing unit normal v can be written as 

i.^^) =±|(-Vv.(z) lf[\~Vv{z) if. (7) 

To see this, note that for x £ dfl the direction of ^{x) is uniquely described by 
the property that for every smooth curve f in dft satisfying ^{0) = x = ( tp(z) ) 



the vectors ^{x) and ^'(0) are orthogonal. Since by ([5]), locally ^(i) 
where C(0) — z, ^ follows from the identity 



vim) 



-v^(^) 1) (v^(c[o))^. c'(0)) ^ ■ + ■ = 
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Proof of Proposition^^ Let x G Oil be arbitrary. Working locally near x, for sim- 
plicity we may assume without loss of generality that there exists (p be as in the 
previous remark, i.e., that no rotation is needed for to be the subgraph of a 
smooth function. Then a: = ( ^fz) ) for some z. Using ([7]), the derivative of T* is 

Here, c = ± | ( ^^vi^) j | ^ ^ 0. In particular, we obtain 

det r*'(z, 0) = c . det L ' , -^^(^r) =c-det(' ^ ^ 

by applying the Gauss-Jordan elimination algorithm. The inverse function theorem 
asserts that T* and hence T is locally a smooth diffeomorphism. Because x e dft 
was arbitrary, all that remains to show is that T is injective if (5 > is small enough. 

By the above argument for every x € dQ there exists an open neighborhood Ox of 
X in dfl and Sx > such that T is a smooth diffeomorphism from Ox x {—dz, Sz) to a 
neighborhood of x in . By compactness of dfl we can choose finitely many such 
Oxi (* = 1, ■ • • , w) which already cover dfl. It is easily proved by contradiction that 
we can find S > such that for every x £ dQ there exists an index k{x) S {1, . . . , m} 
with the property that B4s{x) n dft C Ox^^^^, where Br{a) denotes the open ball 
with center a and radius r. We pick 6 such that S < dx^ for all i = 1, . . . , m. 

For this choice of 6, T is injective. To see this, let T{yi,ti) = T{y2,t2) where 
2/1)2/2 G di} and ii,<2 G {—^7^)- We estimate 

= |T(yi, ii) ~ T{y2, h)\ > \yi - y2\ - {\tMyi)\ + \t2Hy2)\) > I2/1 - 2/2! - 26. 

This shows \yi - 2/2! < 25 and thus 2/2 e -S45(yi), hence yi,y2 G 0^^, fc = fc(2/2)- 
By construction, T is injective on Ok x (—(5,(5), hence yi = ?/2 and = /;2, proving 
the claim. □ 

Lemma 3. The set D := D(A/j) n C°°(ri) is an operator core for Ar, i.e., D is 
dense in Z?(Ar) with respect to the graph norm. 

Proof. As Ar is a generator, the space PlneiN ^(^fl) ^ core for Ar [H Proposi- 
tion II. 1.8]. Moreover, 

i?(i, Ar) n 0(0)) c 

for every m G ff^o by the regularization properties of elliptic operators [51 Re- 
mark 2.5.1.2]. By a standard Sobolev embedding theorem [5, Section 5.6], 

£'(A^) = i?(l, A;^J"C(f7) c H'^'^n) c C2"^[^]"^(fl) 
for all n > ^ . Letting n — > 00 we obtain the assertion. □ 

We remark that for u e C°°(il) the normal derivative exists in the classical sense. 
For these functions, u E D{Aj^) is equivalent to ([2]), and we will frequently use the 
boundary condition in this way. 

Let G2{t) denote the Co-semigroup on L^(R^) with generator 

D{A2) := {u e L2(1R,^)|Au e L\R^)} , 
A2U := Alt. 
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The semigroup G2{t) leaves the space Co(JIl^)ni^(R^) invariant, and its restriction 
extends continuously to a positive, contractive Co-semigroup on Co(Il^), denoted 
by Go{t). The generator of this semigroup is 

D{Ao) := {u e Co(Il^)|A^/ € Co(l!l^)} , 
Aqu :— Au. 

We will refer to both semigroups as the Gaussian semigroup. For more details about 
the Gaussian semigroup, we refer to [U Chapter 3.7]. 



3. Extension Operator 

Given a smooth bounded open set and a smooth function /?: dft 1R+, we 
construct an extension operator Ep which satisfies the assumptions under which 
we will prove ([4]) in Section |4l For (3 — 0, the operator is similar to, but slightly 
simpler than the extension operator in [5^ Section II. 5. 4]. However, the properties 
which we prove here may also be of independent interest. 

For the whole section, let S and T be as in Proposition [TJ We start by fixing a 
"kinking function" g. First choose a function gi having the following properties. 

(a) gi e C°°([0,oo) x [0,oo)) 

(b) < gi(7, t) < 1 for all 7, t > 

(c) gi{j, t) ^0 for alH > f and 7 > 

(d) £'1(7,0) = 1 for aU 7 > 

(e) f£^i(7,0) = -27forall7>0 

(f) ^^?l(7,0) = 472forall7>0 

Here ■§[gi denotes the partial derivative of gi with respect to the second argument. 
For example, we may choose £'1(7, t) := exp {~2^t) x{t) where x is a smooth cut-off 
function such that x = 1 near 0. 
Now define g: fl'-' — > IR to be 



g{x) 



gi{(3{z),t), ifx = T{z,t),0<t<S, 
0, otherwise. 



Note that g is well-defined since T is injective, and it is smooth by construction. 

Definition 4 (Reflection at the Boundary). Let x G T{dn x (~6,S)), x — T{z,t). 
We call X :— T{z,~t) the (orthogonal) reflection of x at the boundary dQ. For a 
function u: ^ R we define the reflected function 

U-.Tidflx {0,5))-^R, u{x):=u{i). 

We define the extension operator Ep belonging to f3 as 

i?^: C(n)^Co(lR^), Epu:^h % (9) 

I gu, on il . 

Here gu is understood to be outside T{dil x (0,(5)) because g equals in that 
region. 

Lemma 5. The operator Ep is well-defined, linear, positive, contractive and an 
extension operator, i.e., REp = /, where R: Co(IIl^) GiVl), u i-^- u\ii. 
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Proof. Let u G C(ri). By property ([d|, the function Epu is continuous on R^. 
Since it has compact support, Epu S Co(R^). Positivity and contractivity foUow 
from property (|b|. The other two properties are obvious from □ 



We now turn towards a more interesting property of Ep: we prove that it maps 
D as defined Lemma [3] into D{Aq). This extensive calculation is split into several 
lemmata. Most calculations will be carried out in local coordinates, i.e., locally at 
xo = T*{zo, 0) € dfl, where we represent all functions with respect to the charts as 
follows. Here T* is defined as in ([6]). 

u*{z,t) := u{T*{z,t)) u*{z,t) := u{T*{z,t)) 

P*{z) :=/3(T*(z,0)) g*{z,t) -.^ g{T*{z,t)) ^ g^{(3* {z),t) 



In the following we will adhere to the usual notation for normal derivatives, i.e., 

I 

du 



^ denotes the directional derivative of g along the outwards pointing unit normal 



with respect to the domain of g. Note that for functions defined on H,'-^ this means 
that ^ — —Vg ■ V, where v always denotes the outwards pointing unit normal of 
VI. 

Lemma 6. Let u e D. Then E^u e £'(A2) and [A2{Eu))\n = Aru. 

Proof. The continuous function Ejsu has compact support, hence EjjU g L^(]R^). 
Moreover, E/^u is smooth away from dfl being the composition of smooth functions. 
Thus the measurable function 

^ _ (Au, on n, 
\A{gu), on Q . 

is defined outside dfl which is a set of measure zero. As u and gu are smooth up 
to dil, f is bounded. Note that / has compact support, hence / £ L^(]R^). 

For the assertion of the lemma, it remains to show that / = A{Epu) in the sense 
of distributions. For this we calculate the (classical) normal derivative of gu using 
that u satisfies For z £ dil we have 

du , . u(z + hi^(z)) — u(z) u{z — hv{z)) — u{z) ^, , , 

— (z) = — lim = — hm = — )(z)u[z) 

dv h^o h h^o h 

and 

^(z) = - lim gl^ + ^'^l^)) - g(^) ^ _ lij^ gi(/3(z),fe) - gi{p{z),0) ^ 
dv h^a h h^a h 

This implies 
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Now let (/? e I?(]R^) be an arbitrary test function. From the above calculations, 
the classical Green formula 21 Section II. 1.3] and Qu\dit — we obtain 



/ [Epu)/S.ip~ / uA(^+ / gu/S.ip 



. f f dif du \ , 

Auip + / u- —if da 

n Jian)+ \ oiy dv J 

A / ~\ f { -d^ d(gu) , , 

A[gu)ip + / gu- — ip da 

nc J {an)- \ oiy dv ' 



Auip+ A{gu) ip^ f(p, 
a Jnc jRN 

where (951)+ is understood as the (oriented) boundary of fi, whereas {d^l)~ de- 
notes the (oriented) boundary of ft'-' . This shows A{Ei3u) = f in the sense of 
distributions. □ 



Remark 7. The above lemma tells us that A(i?^u) is a function. To sec that 
-E/3M e D{Ao), it remains to show that A(£'^u) e Co(]R.^). We already know that 
A{Epu) has compact support and is smooth on \ 957. Thus it suffices to show 
that the function can continuously be extended to d^l. This is a local property. In 
fact, since we already know that the limits from the inside and the outside both 
exist, it suffices to show that Au(a:o) = A{gu){xQ) for every xq € dQ. 

Let Xo £ dfl be fixed. To simplify notation, we may assume that v{xo) = cn 
without loss of generality, exploiting the rotational invariance of the Laplacian. 
Here and in the following, e„ denotes the the n^^ unit vector in R^. Moreover, 
since we treat the problem locally, we may work in local coordinates, a:o — T*{zq, 0). 
We start by calculating the partial derivatives of T*~^, where T* is defined as in ([6]). 

Lemma 8. For n e {1, . . . , N}, 




if n = N. 



Proof. The assumption v{zo) — cn implies V(^(zo) = due to As in ([8]), this 
shows T*'(zo,0) = /. By the inverse function theorem. 

For the partial derivatives at xq this means 
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To calculate the second derivatives, we employ a differentiation rule for matrices, 
i,{A(t)-^) = -A-'(t)A'(t)A-Ht}. 

^-(r"(r-(.)))-'(A(r-(T.-.,.))))(r-(T.-.w))-'.„ 

If we denote the entries of T*' by tij {i, j ~ 1, . . . , N), we can proceed as follows. 

^(t«(T-'w)).vt„(r.-w).(^(T-'))w 

For X ^ xo this yields 



dxn 



(t* ^{xo) \ = Viy(zo,0)e„ = -^Uj{zo,0), 



where for notational simplicity we use as an alias for the variable t. Inserting 
this expression into the above identity, we arrive at 



i,j = l,...,N 



,2 



. Unizo,0) =-j-^T*{zo,0). 

In combination with formula ^ this finishes the proof. □ 

Having the derivatives of the charts at hand, we are able to calculate all deriva- 
tives in local coordinates. 

Lemma 9. Let f and f* be functions such that locally f*{z,t) = f{T(z,t)). Then 
V/(a;o) = Vr(^o,0), 

n=l " n=l " 

/n particular, 

a2 

Vg(xo) = (0 -2/3(a;o)) , ^g{xo) = 4/3(xo)2 + 2/3(a;o) ^ g;j^(^o). 



ri=l 



Proo/. Differentiating f{x) = /* (r*"^(x)) we obtain 



& 



+ (V/-) (t-'w) (t-)) w, 



where Hf* — ( — f* I denotes the Hessian matrix of f*. By using 

Lemma |S] and summing up, we arrive at the desired formulae for x — xq. 
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Concerning g we remark that g*{z,0) — g{j3{z),0) — 1 implies ■^g*{zo,0) — 
{n — 1, . . . , N — 1). On the other hand, the derivatives with respect to t equal 

^e*(^o,0) = ^f?i(/3*(^o),0) = -2/3*(zo) = -m^o), 

— e*(zo,0) = ^^>i(/3*(zo),0) =4/3*(zo)' = 4/3(xo)2. 

With this information, the formulae for g follow from the general formulae. □ 

Finally, it is easy to calculate the relation between the derivatives of the function 
and its reflection at the boundary in local coordinates. It suffices to observe that 

t) — u{T[z, t)) = u{T{z, —t)) — u* {z, — i). 

From this we deduce the following formulae. 

u* (z, t) ~ u* {z, —t) 

d d d d 

——u*{z,t) = ——u*{z,-t) —u*{z,t) = -—u*{z,-t) 

azn ozn at ot 

Q2 Q2 Q1 Q2 

Now we are ready to prove continuity of l^[Epu) at xq. 

Proposition 10. For every u ^ D, Au{xo) — A{gu){xo) . 

Proof. Note that 

Ou 

—u*{za, 0) = -—u*{zq, 0) = -^(a;o) = P{xo)u{xo) = f3{xo)u{xo). 
ot ot Ov 

We use the formulae of this section to obtain the desired identity. 
A.(gu){xQ) = A.g{x(3) u{xq) + 2'Vg{xo) ■ Vuixa) + g{xQ) Au{xq) 

^-1 Q2 Q 

= AP{xofn{xQ) + 2l3{xo)u{xo) ^ ^'/'(^o) - 4:P{xo)-^u* (zq, 0) 

n=l " 

Q2 Q2 Q N^l 

+ ^"*(^o,0) + ^r(zo,0)- ^u*(zo,0) ^"^(^o) 

n=l " n=l " 

N-l q2 q2 g N~l ^2 

n— 1 ^ n— 1 ^ 



AM(a;o)- 



□ 



The following theorem is the main result of this section. As explained in Re- 
mark [71 it follows by combining Lemma [6] and the last proposition. Even though 
Theorem 111 I is also true for the usual extension operator for Lipschitz domains [10| 
VI. §3, Theorem 5], that operator fails to be contractive and thus is more difficult 
to handle for the application in Section [4] 

Theorem 11. The operator Ep maps D into _D(Ao). 

Corollary 12. The operator Ep maps D{Aji) into D{Ao). 
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Proof. There exists a constant C > satisfying ||-E/3u||£)(Ao) ^ C||'"I|d(Ab) for all 
u £ D. To see this, note that on 

\\A{ug)\\^ = \\Aug + 2VuVg + uAg\\^ 

< IIAulU ||^?IL + 2 IIV^^II^ (e IIA^IU + II^IL) + II^IL II A^'IL 

for every e > 0. Similarly, < and ||A{t||^ < C (||Am||^ + using 

the definition of {t as a composition of u and a function involving T, where C > 
depends only depends on a estimate on the derivatives of T. Noting that g and its 
derivatives are bounded by assumption, we sec that there exists a C > as claimed. 

As 13 is a core of Ar, the above estimate shows that there exists a unique 
continuous extension of -E/jId to D(Afi), and that this operator still takes values in 
£'(Ao). Because D{Afj) is continuously embedded into C{fl) and Ep is continuous 
on C(r2), this extension agrees with -E/3|£i(An)- Thus the claim is proved. □ 

4. Approximation Result 

In this section, we we prove that if Ep: C(S1) Co(R^) is a contractive exten- 
sion operator mapping an operator core D for Aj^ into D{Aq), then formula ^ 
holds. Note that the operator defined in ([9]) has this properties as shown in the 
preceding section. The tool we use for the proof is the following approximation 
result for semigroups due to Chernoff. 

Theorem 13 C^, Theorem in.5. 2]). Let X be a Banach space. Consider a function 

V : [0,oo) ^^(A:) 

satisfying V{0) = / and \\V{t)"'\\ < M for all t > 0, m e K and some M > 1. 
Assume that 

Vih)x-x 
Ax hm 

h^O h 

exists for all on x Cz D d X , where D and (Ag — A)D are dense subspaces in X for 
some Aq > 0. 

Then {A, D) is closable and A generates a bounded Go-semigroup T{t), which is 
given by 

T{t)x^ lim(n^))"x 

for every x Cz X locally uniformly with respect to t > 0. 

We apply the theorem by setting 

X := Cp, Vit) := RGo{t)Ep, A A^. (10) 

As D is an operator core for A/j, the density conditions are fulfilled because A — Ar 
is an isomorphism between D{Aii) with the graph norm and C(ri) for every A > 0. 

Theorem 14. Let Ep: C(f2) — > Co(JIl^) be a contractive extension operator which 
maps an operator core D for Ar into I3(Ao). Then formula (jj]) holds true. 

Proof. We check the conditions of Chernoff's product formula with the choices 
made in pU)) . The fact 1^(0) = / is equivalent to Ep being an extension operator. 
Since all three of their factors are contractions, the operators V{t) are contractions 
for every i > 0, thus < 1; in particular V{t) is a bounded operator for 

every t > Q. The density assumptions on D are fulfilled because D is an operator 
core for Ar. 
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Now let u E D be arbitrary. By assumption, Epu G D{Ao). By definition of the 
infinitesimal generator, 

h h 

Since the function Epu agrees with u on fl, they represent the same distribution 
acting on the test functions I?(r2). This means that they have the same distribu- 
tional derivatives, hence RAoE/ju = Arm. Having checked all the conditions of 
Theorem [T3l we deduce that indeed (|4]) holds true. □ 

Corollary 15. Formula Q holds true for the operator Ep defined in ([9|. 

Remark 16. As a special case, we may choose /3 = 0. Then Ar — Ajv is the Lapla- 
cian with Neumann boundary conditions. In this case, Eq is the reflection without 
"kinking" , corresponding to certain numeric schemes where Neumann boundary 
conditions are realized as in ([3]) . A different extension operator for Neumann bound- 
ary conditions would be given by extending constantly along the outwards pointing 
unit normal and again multiplying by a cut-off function. This corresponds to a 
first-order accurate boundary condition approximation, see again |1H Section 8.3]. 
Although this might seem more natural at first, it is not obvious whether formula @ 
is true for this choice of i?^. Unfortunately, Chernoff's theorem cannot be applied 
again because AEfju fails to be continuous at dfl as can easily be seen. 

5. DiRiCHLET Boundary Conditions 

Next we treat the model problem of an elliptic operator on a bounded set, the 
Laplacian with Dirichlct boundary conditions. Typically, all results about ellip- 
tic operators are much simpler for this special case. Surprisingly, for the aim of 
this article there arise completely different problems than for Robin and Neumann 
boundary conditions. This is the reason why we consider it worthwhile to treat this 
operator in detail. 

The Laplacian with Dirichlct boundary conditions defined by 

D{Ad) := {u e Co(f^)|Au e Co{n)} , Adu := Au 

generates an positive, contractive Co-semigroup To{t) on Co(ri) [P, Theorem 6.1.8]. 
Formally, the boundary conditions ([2]) become the Dirichlct boundary conditions 
M =: on dVl in the limit (3 oo. This observation can be made precise, cf. [T31 
Proposition 3.5.3]. As we want to prove an analogue of ([4]) for To{t), we have to 
define an appropriate extension operator Eao for P = oo. Taking the limit in ([9]), 
we arrive at 

Eoc - Co(0) ^ Co(]R^), E^u:= 

Note that we had to replace C{n) by Co(l^) as we require EooU to be continuous. 
Unfortunately, we cannot simply replace Efs by E^o in formula ^ because the 
iteration scheme does not remain in Co(f^), hence leaving the domain of i?oo- 

However, the analogue formula is well-defined (and true) in the L^-context. To 
see this, note that L^(ri) is a closed subspace of L^(R^) if we consider its functions 
to be extended by zero. Then the identity mapping takes the role of Eao, and the 
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restriction becomes multiplication with Iji- Thus, the analogue of formula (jlj for 
Dirichlet boundary conditions reads 

TD.2{th = lim (laG2(:^))"w in L^{n) for every u £ L^{n), (11) 

where denotes the Dirichlet semigroup on L^(ri) generated by the Laplacian 
on L^(ri) with domain i?o(ri) n H'^{n). Indeed, formula (fTT|) remains true even if 
fl has merely Lipschitz regular boundary, cf. [8]. 

It is interesting to note that (|11|) cannot be proved using Chernoff's product 
formula in the way we did in SectionlH For this, a dense subspace of Hq (fi) fliJ^ (fi) 
would have to be contained in H^{R^), where both spaces carry the graph norm of 
the Laplacian. But then, continuity asserts i?o(ri) n H^{il) C H^(R^). However, 
this cannot be true. In fact, a function in Co{n) D C°°{i}) C Hq{Q) n H'^{Q) whose 
normal derivative does not vanish is not an element of iJ^(]R^). 

Despite those problems, it is possible to prove a similar result in the same spirit 
as in Section |4] even in Co(r2). For this, we need to replace In by a sequence 
of smooth interior cut-off functions. But we have to assure that they exhaust f2 
sufficiently fast compared to the decay of functions in a core for A^. So we start 
by a investigation of that decay. 

Lemma 17. Given a Dirichlet regular bounded set 51, there exists m Cz K having 
the following property. Given t > 0, there exists a neighborhood Ut of dil such that 
the estimate 

\u{x)\<t^Ul-Anru\\^ 
holds for every x Cz Q HUt and every u G D(N^). 

Proof. It is well-known that To{t) has a kernel representation with a continuous 
non-negative symmetric kernel ks{x,y) which vanishes on dfl and is dominated by 
the Gaussian kernel [3l Section 3.4]. Let m > -j be fixed. The integral formula for 
powers of the resolvent f4l Corollary 2.1.11] shows that (/ — Ad)"™ is a positive 
kernel operator with the continuous non-negative symmetric kernel 

poo 1 

Jo [m-iy. 

which vanishes for x G dfl. Using compactness of fl and d^l we deduce that for any 
e > there exists a neighborhood of 9f2 such that x € implies k{x, y) < e for 
all y G ri. Define Ut ■= S^, where e := -jl^i-. 

Now fix u G £'(AS) and define v:= {I - Ad)'"m G Co(ll). For x eQnUt, i.e., 
a; G 51 n S'e, we obtain 



\u{x)\ <{I- Az,)-™ \v\ ^ / kix,yMy)\ dy < e \n\ \\v\\^ = t' - A^)'»u|U ■ 

Jn 

This concludes the proof. □ 



We have already explained why we cannot use i?oo as extension operator. In- 
stead, we choose 

Ed: Co(fl) -> CoCR"^), Edu := 
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similarly to Here g denotes a cut-off function that equals 1 near d^l. Using 
the same ideas as in Section [3] it can be shown that Ed is a contractive extension 
operator that maps Z?(Ad) n C°°(il) into D{Ao). In fact, the main difference to 
Section[3]is that we know Au e Co(ri) for u e D{Ad) which makes it easy to check 
the continuity of A(E]ju), significantly shortening the chain of arguments. 

Now let m be as in the above lemma, and choose a family {Ut)t>o as in the 
lemma. For every t > we fix a suitable cut-off function xt G Co(r2) satisfying 
< xt < 1; and Xt{x) — 1 if X e n\Ut. Moreover, define xo ■— 1^- simplify 
notation, we use the multiplication operator xt as an operator from Co(IEl^) to 
Co(0) and xo as the restriction from Co(R^) to C(ri), whenever they are applied 
to functions in Co(Il^). 

We remark that in view of the kernel of Tjj (t) {t > 0) being strictly positive in 
the interior of due to the strong maximum principle, it can be seen that for every 
compact set K C ft there exists to > such that Ut and K are disjoint whenever 
t < io, implying that Xt ^ lj2 pointwise as i — > 0. In this sense, the next result is 
another flavor of formula (fTTj) . 

Theorem 18. Let m G IN and {xt)t>a be as above. Then 

TD{t)u^ lim {x±Go{^)EDyu 

for every u G Co(ri) uniformly on [0,T] for every T > 0. 

Proof. We apply Theorem [T3l to the operators 

V{t) : Com ^ Co(r!), u ^ xtGo{t)EDU. 

The properties V^(0) = / and < 1 for every t > Q and n G IN are obvious 

from the properties of Ed and the Gaussian semigroup. Let D := D{A!^)r]G°° (fl), 
which is a core for Aq. This choice makes the density conditions automatic once 
we show that the limit operator is Aq. 

It only remains to prove the convergence to Au on D. For this, let u G -D. In 
particular u G ^'(Ad), thus Au G Co(f^). Note that 



V{t)u- 



Au 



XtGoit)EDU 



< 



Xt 



Go{t)EDU-EDU 



t 



t 

Au 



du — u 


- Am 








OO 




+ 


XtEnu - u 


) 


OO 


t 



XtAu - Au| 



We estimate the three summands separately. After estimating xt by 1 in the 
first expression, convergence to zero follows from E^u G Z?(Ao) and the fact that 
Ao{Edu) = Au on 17. The third summand can be estimated by sup^^jj^ 2 |Au(a::)| 
using that Xt = 1 on \ J7t . But since we assumed that Ut leaves any compact set 
C O for small t, this expression becomes small as t because Au G Co(ri). 
The second summand can be estimated with the help of Lemma [T71 We obtain 



XtEnu 



t 



< - sup |u(a;)| < 2t\\{I- Aof 
t xeUt 







as t ^ 0. Together, these three estimates show the convergence of the difference 
quotient to Au as t tends to zero. We have checked the assumptions of Chernoff's 
product formula, thus proving the claim of the theorem. □ 
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6. Conclusion 

It is a direct consequence of (jlj that Tii{t) is a positive semigroup. Because 
the operators on the right are L°°-contractive, it is also clear that Tfi{t) is L°°- 
contractive, thus submarkovian. In the same way other properties of the limiting 
semigroup can be deduced by such an approximation formula, as long as they are 
preserved when taking limits in the strong operator topology. To obtain further 
properties, it might help to modify the formula a little bit. 

So far, we have only considered the Gaussian semigroup as underlying tool. 
However, it can be seen from the proofs that actually we used only few properties of 
the Gaussian semigroup. More precisely, we only used that Go(t) is a contraction on 
Co(IIl^) and that any continuous function u such that the support of u is contained 
in a given neighborhood of il and Am is continuous on is in the domain of 
the generator of Go{t). Thus we could replace Go{t) with other semigroups, for 
example with the semigroup generated by the Laplacian with Dirichlet or Neumann 
boundary conditions, on a larger bounded set fi' C . Then (|4]) becomes an 
approximation formula where the approximating operators are compact. Note, 
however, that this does not imply that Tfi{t) is compact, as the limit is only in the 
strong operator topology. 

Similarly, we can try to approximate Tjiit) only in terms of operators on C(f7), 
i.e., without any extension to R^, to obtain an intrinsic approximation. The most 
natural Trotter-like candidate of this kind would be 5„(i) := {Td{^)Tn{^))" , 
where a and /3 are positive numbers such that a + f3 — 1 and TD{t) and T^lt) 
denote the Dirichlet and Neumann semigroups on C(f2), respectively. It is known, 
however, that lim S'„(i) = T]j(t) in the strong operator topology on L^(r2), whenever 
a > 0, see [7]. 

Recall that regarding the extension operator we used only two of its properties 
in Section |4l namely contractivity and some regularity of the extended function. 
By definition of the extension operator, contractivity came for free. This is due to 
the rather special definition of Efj and the choice of spaces and is a very convenient 
prerequisite for the application of Theorem 1131 although not a necessary one. 

Assume that we replace Ep by some other, non-contractive extension operator 
E. This is a natural consideration because most extension operators are non- 
contractive. In fact, it is easy to see that no operator extending C^(ri) to C^(R^) 
can be contractive. This also shows that it is a very special property for an extension 
operator to be contractive and to preserve the regularity of functions in D(Ar). 

For such an extension operator E^ it is considerably more difficult to check 
whether [RG{^)E^ is uniformly bounded in operator norm with respect to n. 
Because it is hard to control such iterated applications of the Gaussian semi- 
group, one could try to estimate each factor separately. Then one has to show 
that ||i?G'o(i)£'|| <l + ct for some c > 0, leading to the upper bound e'^*. The short 
time diffusion through the boundary, however, is of order 0{-\/t), see [9]. This is 
why in general only estimates of the kind \\RGQ[t)E\\ = l + 0{^/t) can be obtained. 

Almost the same reasoning applies if C(0) is replaced by an L^'-space, for exam- 
ple by Z/^(f2). As is bounded, uniform convergence already implies convergence 
in L^(ri), hence 

TB{t)u = lim {RG^{^)EaYu in L^{n) for every u G CfH) 
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by what we have aheady shown. If we want to extend this result to u e L^(f2), it 
suffices to show that the approximating operators (i?Go(^)i?/3) remain bounded in 
the norm of operators on L^{fl). Here again, there arise difficulties which are similar 
to those mentioned in the preceding paragraph because no non-trivial extension 
operator from L^(r2) to L^(R^) is contractive. This shows that for our applications 
the space C{^1) has significant advantages. 

A related question is whether (j4]) remains true if the assumption /3 > is 
dropped. We mention that it can be seen that for any /3 e L°°{dil) the Lapla- 
cian with Robin boundary conditions is the generator of a semigroup on (il) and 
thus this question makes sense. But if we define Efj as in Section [31 we do not 
even in C(n) obtain a contraction if < for a point z E dil, causing the same 
problems again. Moreover, it is clear that the assumption || (i?Go(i)£';3)"|| which is 
needed for Theorem[T3] cannot be fulfilled since the candidate limit semigroup Tii{t) 
will not be bounded. But the latter is merely a problem of rescaling, compare [H 
Corollary III.5.3]. 

It should be possible to extend the results to smooth unbounded open sets il 
without difficulties because most arguments are local. However, the other calcu- 
lations become even more technical. This is why we have restricted ourselves to 
bounded domains. 

On the other hand, choosing a different (contractive) extension operator will usu- 
ally change the situation completely. For example. Theorem [13] cannot be applied 
for the constant extension as in Remark [TH reflecting the fact that a worse numer- 
ical approximation of the normal derivative leads to worse convergence behavior. 
But that extension operator can be defined even for convex domains without any 
smoothness assumptions, which might provide an alternative approximation scheme 
for less smooth domains. It is easy to come up with various other extension opera- 
tors when trying to find an approximation formula such as ^ for (not necessarily 
convex) sets with non-smooth boundary. This is ongoing work and might be the 
topic of a future publication. 

References 

1. W. Arendt, C. Batty, M. Hicbcr, and F. Ncubrander, Vector- Valued Laplace Transforms and 
Cauchy Problems, Birkhauser, 2001. 

2. R. Dautray and J.-L. Lions, Mathematical Analysis and Numerical Methods for Science and 
Technology 1: Physical Origins and Classical Methods, Springer- Verlag Berlin, 1990. 

3. E.B. Davies, Heat kernels and spectral theory, vol. 92, Cambridge Tracts in Mathematics, 
1989. 

4. K.-J. Engel and R. Nagel, One-Parameter Semigroups for Linear Evolution Equations, 
Springer, 2000. 

5. L.C. Evans, Partial Differential Equations, American Mathematical Society, 1998. 

6. P. Grisvard, Elliptic Problems in Nonsmooth Domains, Pitman, Boston, 1985. 

7. T. Kato, Trotter's product formula for an arbitrary pair of self-adjoint contraction semigroups 
Topics in functional analysis. Adv. Math. Suppl. Stud 3 (1978), 185-195. 

8. M. Matolcsi and R. Shvidkoy, Trotter's product formula for projections, Archiv der Mathe- 
matik 81 (2003), no. 3, 309-317. 

9. M. Miranda Jr, D. Pallara, F. Paronetto, and M. Preunkert, Short-time Heat Flow and Func- 
tions of Bounded Variation in K,^, Annales, faculte des sciences toulouse mathematiques 16 
(2007), no. 1, 125. 

10. E.M. Stein, Singular Integrals and Differentiability Properties of Functions, Princeton Uni- 
versity Press, 1970. 



16 



R. NITTKA 



11. J.C. Strikwerda, Finite Difference Sehernes and Partial Differential Equations, Society for 
Industrial Mathematics, 2004. 

12. H.F. Trotter, On the Product of Semi-Groups of Operators, Proceedings of the American 
Mathematical Society 10 (1959), no. 4, 545-551. 

13. M. Warma, The Laplacian with General Robin Boundary Conditions, Ph.D. thesis, University 
of Ulm, 2002. 

14. , The Robin and Wentzell-Robin Laplacians on Lipschitz Domains, Semigroup Forum 

73 (2006), no. 1, 10-30. 

University of Ulm, Institute of Applied Analysis, Germany 
E-mail address: robin.nittka9uni-ulm.de 



